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How to approach continuum physics in lattice Weinberg - Salam model 
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We investigate lattice Weinberg - Salam model without fermions numerically for the realistic 
choice of coupling constants correspondent to the value of the Weinberg angle 8w ~ 30°, and 
bare fine structure constant around a ~ y^. We consider the values of the scalar self coupling 
corresponding to Higgs mass Mh ~ 100, 150, 270 GeV. It has been found that nonperturbative 
effects become important while approaching continuum physics within the lattice model. When the 
ultraviolet cutoff A = — (where a is the lattice spacing) is increased and achieves the value around 1 
TeV one encounters the fluctuational region (on the phase diagram of the lattice model), where the 
fluctuations of the scalar field become strong. The classical Nambu monopole can be considered as 
an embryo of the unphysical symmetric phase within the physical phase. In the fluctuational region 
q ■ quantum Nambu monopoles are dense and, therefore, the use of the perturbation expansion around 

^s^j ' trivial vacuum in this region is limited. Further increase of the cutoff is accompanied by a transition 

to the region of the phase diagram, where the scalar field is not condensed (this happens at the 
^ , value of A around 1.4 TeV for the considered lattice sizes). Within this region further increase of 

the cutoff is possible although we do not observe this in details due to the strong fluctuations of the 
gauge boson correlator. Both mentioned above regions look unphysical. Therefore we come to the 
conclusion that the maximal value of the cutoff admitted within lattice Electroweak theory cannot 
exceed the value of the order of 1 TeV. 



PACS numbers: 12.15.-y, 11. 15. Ha, 12.10.Dm 

in: 

i 1, I. INTRODUCTION 

It is well - known 1] , that the finite temperature perturbation expansion breaks down at the temperatures above the 
\Q ' electroweak transition/crossover already for Higgs masses above about 60 GeV. Therefore the present lower bound 
on the Higgs mass requires the use of nonperturbative techniques while investigating electroweak physics at high 
" temperature. 

Nambu monopoles are not described by means of a perturbation expansion around the trivial vacuum background. 
Therefore, nonperturbative methods should be used in order to investigate their physics. However, their mass is 
C> estimated at the Tev scale. That's why at zero temperature and at the energies much less than 1 Tev their effect on 
physical observables is negligible. However, when energy of the processes approaches 1 Tev we expect these objects 
influence the dynamics. Recently the indications in favor of this point of view were indeed found 

m. 

In this paper we consider lattice realization of zero temperature Electroweak theory (without fermions). The phase 
diagram of the correspondent lattice model contains physical Higgs phase, where scalar field is condensed and gauge 
bosons Z and W acquire their masses. This physical phase is bounded by the phase transition surface. Crossing this 
surface one leaves the Higgs phase and enters the phase of the lattice theory, where the scalar field is not condensed. 

In the lattice theory the ultraviolet cutoff is finite and is equal to the momentum A = — (see, for example, |5|), 
where a is the lattice spacing. The physical scale can be fixed, for example, using the value of the Z-boson mass 
A/P hys _ 90 GeV. Therefore the lattice spacing is evaluated to be a ~ [90 GeV] l Mz, where Mz is the Z boson 
mass in lattice units. Within the physical phase of the theory the lines of constant physics (LCP) are defined that 
correspond to constant renormalized physical couplings (the fine structure constant a, the Weinberg angle 9yv, and 
Higgs mass to Z-boson mass ratio rj = Mh/Mz). The points on LCP are parametrized by the lattice spacing. Our 
observation is that the LCP corresponding to realistic values of a, 9w, and rj crosses the transition between the two 
"phases" at a certain value a = a c and for a < a c the scalar field is not condensed. We denote the corresponding 
value of the cutoff A c = — . Our estimate for the considered values of the Higgs mass Mjj ~ 100, 160,270 Gev is 
A c = 1.4 ±0.2 Tev (for the considered lattice sizes). We do not observe the dependence of A c on the lattice size. That's 
why the value A c might appear as the maximal possible value of the cutoff allowed in the conventional Electroweak 
theory. 
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It is important to compare this result with the limitations on the Ultraviolet Cutoff, that come from the perturbation 
theory. From the point of view of perturbation theory the energy scale 1 TeV appears in the Hierarchy problem 
[f|. Namely, the mass parameter /j? for the scalar field receives a quadratically divergent contribution in one loop. 
Therefore, the initial mass parameter (/i 2 = —X c v 2 , where v is the vacuum average of the scalar field) should be set to 
infinity in such a way that the renormalized mass [i 2 R remains negative and finite. This is the content of the so-called 
fine tuning. It is commonly believed that this fine tuning is not natural |6j and, therefore, one should set up the finite 
ultraviolet cutoff A. From the requirement that the one-loop contribution to /s 2 is less than 10|/i^| one derives that 
A ~ 1 TeV. However, strictly speaking, the possibility that the mentioned fine tuning takes place is not excluded. 

In the perturbation theory there is also more solid limitation on the Ultraviolet cutoff. It appears as a consequence 
of the triviality problem, which is related to Landau pole in scalar field self coupling A and in the fine structure 
constant a. The Landau pole in fine structure constant is related to the fermion loops and, therefore, has no direct 
connection with our lattice result (we neglect dynamical fermions in our consideration) . Due to the Landau pole the 
renormalized A is zero, and the only way to keep it equal to its measured value is to impose the limitation on the 
cutoff. That's why the Electroweak theory is usually thought of as a finite cutoff theory. For small Higgs masses (less 
than about 350 GeV) the correspondent energy scale A" calculated within the perturbation theory is much larger, 
than 1 Tev. In particular, for Mh ~ 300 GeV we have A° ~ 1000 TeV. It is worth mentioning that for A — > oo the 
perturbation expansion in A cannot be used. In this case Higgs mass approaches its absolute upper bound 35], and 
both triviality and Hierarchy scales approach each other. 

From the previous research we know that the phase diagram in the j3 - 7 plane of the lattice SU(2) Gauge - Higgs 
for any fixed A resembles the phase diagram for the lattice Weinberg - Salam model. The only difference is that in the 
SU(2) Gauge - Higgs model the confinement-deconfinement phase transition corresponding to the U(l) constituents 
of the model is absent. The direct measurement of the renormalized coupling /3r shows [7H20j| that the line of constant 
renormalized coupling constant (with the value close to the experimental one) intersects the phase transition line. 
Also we know from the direct measurements of Mw in the SU(2) Gauge - Higgs model that the ultraviolet cutoff is 
increased when one is moving along this line from the physical Higgs phase to the symmetric phase. 

On the tree level the gauge boson mass in lattice units vanishes on the transition surface at small enough A. This 
means that the tree level estimate predicts the appearance of an infinite ultraviolet cutoff at the transition point 
for small A. At infinite A the tree level estimate gives nonzero values of lattice masses at the transition point. Our 
numerical investigation of SU(2) <Ei U(l) model (at A = 0.0025,0.009,0.001) and previous calculations in the SU(2) 
Gauge Higgs model (both at finite A and at A = 00) showed that for the considered lattice sizes renormalized masses 
do not vanish and the transition is either of the first order or a crossover. (Actually, the situation, when the cutoff 
tends to infinity at the position of the transition point means that there is a second order phase transition.) The 
dependence on the lattice sizes for the SU(2) Gauge Higgs model was investigated, for example, in [l6j]. Namely, for 
P = 8, A ~ 0.00116, where Mh ~ Mw, the correlation lengths were evaluated at the transition points. For different 
lattice sizes (from 12 3 x 28 to 18 3 x 36) no change in correlation length was observed [l6j . 

In table 1 of [2| the data on the ultraviolet cutoff achieved in selected lattice studies of the SU(2) Gauge Higgs 
model are presented. Everywhere (3 is around (3 ~ 8 and the renormalized fine structure constant is around a ~ 1/110. 
This table shows that the maximal value of the cutoff A = - ever achieved in these studies is around 1.4 Tev. 

a 

Thus the predictions on the value of A c given by our lattice study and on the value A° given by the perturbation 
theory contradict with each other. A possible explanation of this contradiction we suggested in Namely, it 
was demonstrated that in the vicinity of the transition there exists the fluctuational region. Within this region the 
application of the perturbation theory is limited. This situation is similar to that of some phenomenological models 
that describe condensed matter systems [36j], where there exists the vicinity of the finite temperature phase transition 
that is also called fluctuational region. In this region the fluctuations of the order parameter become strong. The 
contribution of these fluctuations to certain physical observables becomes larger than the tree level estimate. Thus 
the perturbation theory in these models fails down within the fluctuational region. 

We find that there exists the vicinity of the phase transition between the Higgs phase and the symmetric phase in 
the Weinberg - Salam model, where the fluctuations of the scalar field become strong and the perturbation expansion 
around trivial vacuum cannot be applied. According to the numerical results the continuum theory is to be approached 
within the vicinity of the phase transition, i.e. the cutoff is increased along the line of constant physics when one 
approaches the point of the transition. That's why the conventional prediction on the value of the cutoff admitted in 
the Standard Model based on the perturbation theory may be incorrect. 

In the present paper we proceed the investigation j4j of the model at the value of the scalar self coupling A = 0.009 
(corresponds to the Higgs boson mass around 270 Gev in the vicinity of the phase transition), bare Weinberg angle 
8yy = 30°, and bare fine structure constant around 1/150. The results presented now correspond to essentially larger 
lattices than that of used in Namely, in main results correspond to lattices 8 3 x 16; some results were checked 
on the lattice 12 3 x 16; two points were checked on the lattice 16 4 . Now our main results are obtained on the lattice 
16 4 while the results at the transition point were checked on the lattice 20 3 x 24. 
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In addition we investigate the model at the value of the scalar self coupling A = 0.0025, bare Weinberg angle 
9 W — 30°, and bare fine structure constant around a ~ 1/150. These values of couplings correspond to the Higgs 
boson mass around 150 Gev in the vicinity of the phase transition. The results are obtained using lattices 8 3 x 16, 
12 3 x 16, and 16 4 . We also present results for A = 0.001, 9w = 30°, ao ~ 1/150. These values of couplings correspond 
to the Higgs boson mass around 100 Gev. The results are obtained using lattices 8 3 x 16, 12 3 x 16. 

It is worth mentioning that far from the transition point the renormalized fine structure constant slowly approaches 
the tree level estimate. Contrary to the maximal value of the cutoff the renormalized fine structure constant depends 
on the lattice size. And for the larger lattice the value of a.R is closer to the tree level estimate than for the smaller 
one. For example, for (3 = 12, 7 ~ 1, A = 0.001 (far from the transition point) on the lattice 8 3 x 16 the value of an 
is around 1/130 while on the lattice 12 3 x 16 it is around 1/140. Within the fluctuational region the deviation from 
tree level estimate becomes essentially strong. For example, for A = 0.009,7 = 0.274 (near the transition point) the 
renormalized value of an calculated on the lattice 8 3 x 16 is around 1/99 while on the lattice 20 3 x 24 its value is 
around 1/106. As it is seen from our numerical results and as it will be explained in the Conclusions we guess the 
mentioned finite volume effects present in the value of renormalized a do not affect the main observables we considered 
like the value of A c and the Nambu monopole density. 

We calculate the constraint effective potential V(|$|) for the Higgs field $. In the physical Higgs phase this potential 
has a minimum at a certain nonzero value <f> m of |$|. This shows that the spontaneous breakdown of the Electroweak 
symmetry takes place as it should. However, there exists the vicinity of the phase transition, where the fluctuations of 
the Higgs field are of the order of <f> m while the hight of the "potential barrier" [37| H = V(0) — V(4> m ) is of the order 
of V(<j> m + $4>) ~ Vifim), where 6(j> is the fluctuation of |<&|. We expect that in this region the perturbation expansion 
around trivial vacuum $ = ((f> m ,0) T cannot be applied. This region of the phase diagram is called the fluctuational 
region (FR). 

The nature of the fluctuational region is illustrated by the behavior of quantum Nambu monopoles (U [23| . We 
show that their lattice density increases when the phase transition point is approached. Within the FR these objects 
are so dense that it is not possible at all to speak of them as of single monopoles (38|. Namely, within this region 
the average distance between the Nambu monopoles is of the order of their size. Such complicated configurations 
obviously have nothing to do with the conventional vacuum used in the continuum perturbation theory. 



II. THE LATTICE MODEL UNDER INVESTIGATION 

The lattice Weinberg - Salam Model without fermions contains gauge field U — (U, 9) (where U £ SU(2), e £ 
U(l) are realized as link variables), and the scalar doublet $ Q , (a = 1, 2) defined on sites. 
The action is taken in the form 

S = (3 ((l-|TrL/ p ) + r ^-(l-cos^)) + 

plaquettes ^ ^ 

- 7 ]T Re{<$>+U xy e^<$>) + J2(\^\ 2 + H\^\ 2 ~ I) 2 ), (1) 

xy x 

where the plaquette variables are defined as U p = U xy U yz U^ z U xw , and 9 p = 9 xy + 9 yz — 9 WZ — 9 XW for the plaquette 
composed of the vertices x, y, z, w. Here A is the scalar self coupling, and 7 = 2k, where k corresponds to the constant 
used in the investigations of the SU(2) gauge Higgs model. 9w is the Weinberg angle. 

Bare fine structure constant a is expressed through (3 and 9w a,s a — Kp^A-tf^e-w) • ^ n or der to demonstrate this we 
consider naive continuum limit of ([T]). We set 

JJ x — e iA„(x)a^ gi0 x , M _ e iB ll (x)a ^) 

Here a is the lattice spacing. The field = where - is the conventional U(l) field while is the conventional 
SU(2) field. In continuum limit (fTJ) must become 

J y2 i> 3 yi i>j 

Here Fij = d L Bj — djBi = 2{diBj — djBi) = 2Fij, Gij = diAj — djAi — i[Ai,Aj]. We also have the following 
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correspondence between the plaquette variables and the field strengths: 

TrU x>MV = Tr[l -±Gl v a% 
N 2 

cos N6 X ^ = [1-—F 2 u a 4 ] (4) 

Now in order to clarify the correspondence between constants gi t 2 and /3 we must substitute the expressions for the 
field strengths to (fT]) and compare it to We have: 

i = irir x/3 ' ^ = fj/A (5) 

g l 4tg Vw 92 

Thus 

9i 



tg9 w 

92 



a 



4tt 4tt tt/3(1 + tg 2 9 w ) 



tg 2 e w (6) 



We consider the region of the phase diagram with (3 ~ 12 and 9w ~ 7r/6. Therefore, bare couplings are sin 2 



'w 



0.25; a ~ -j-go ■ These values are to be compared with the experimental 

onessin 2 6W(i00Gev) ~ 0.23; a(100Gev) - 4g. 
The simulations were performed on lattices of sizes 8 3 x 16, 12 3 x 16. For A = 0.0025, 0.009 we investigate the 
system on the lattice 16 4 . The transition point at A = 0.009 was checked using the larger lattice (20 3 x 24). In order 
to simulate the system we used Metropolis algorithm. The acceptance rate is kept around 0.5 via the automatical 
self - tuning of the suggested distribution of the fields. At each step of the suggestion the random value is added 
to the old value of the scalar field while the old value of Gauge field is multiplied by random SU(2) <g> U(l) matrix. 
We use Gaussian distribution both for the random value added to the scalar field and the parameters of the random 
matrix multiplied by the lattice Gauge field. We use two independent parameters for these distributions: one for the 
Gauge fields and another for the scalar field. The program code has been tested for the case of frozen scalar field. 
And the results of the papers Q are reproduced. We also have tested our code for the U(l) field frozen and repeat 
the results of [29j. Far from the transition point the autocorrelation time for the gauge fields is estimated as about 
Na U to ~ 500 Metropolis steps. In the vicinity of the transition point the autocorrelation time is several times larger 
and is about N^ ut0 ~ 1500 Metropolis steps. (The correlation between the values of the gauge field is less than 3% 
for the configurations separated by N^ uto Metropolis steps. Each metropolis step consists of the renewing the fields 
over all the lattice.) The autocorrelation time for the scalar field is essentially smaller than for the gauge fields and 
is of the order of N^ uto ~ 20. The estimated time for preparing the equilibrium starting from the cold start far from 
the phase transition within the Higgs phase is about 18000 Metropolis steps for the considered values of couplings. 
At the same time near the phase transition and within the symmetric phase the estimated time for preparing the 
equilibrium is up to 3 times larger. 



III. THE TREE LEVEL ESTIMATES OF LATTICE QUANTITIES 

At finite A the line of constant renormalized a is not a line of constant physics, because the mass of the Higgs 
boson depends on the position on this line. Thus, in order to investigate the line of constant physics one should vary 
A together with 7 to keep the ratio of lattice masses Mh/Mw constant. 

In order to obtain the tree level estimates let us rewrite the lattice action in an appropriate way. Namely, we define 
the scalar field $ = y^F- We have: 



S = f 3 E ((l-^TrU p ) + -^—(l-co S e p )) + 

plaqucttes 

+ - U xy e ie °»%\ 2 + ^(a^I 2 + A|$ x | 4 ) + w, (7) 

xy x 

where fi 2 = —2(4 + (2A — l)/7), A = 4^-, and to — XV. Here V — L 4 is the lattice volume, and L is the lattice size. 
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For negative u 2 we fix Unitary gauge $2 = 0, Im$i = 0, and introduce the vacuum value of $: v = -^=. We also 
introduce the scalar field cr instead of $1 = v + a. We denote V xy = (U^e l6xy — 1), and obtain: 

5 = /? ]T ((l-iTr[/ p ) + -^— (l-cos6y) + 

plaquettes 

+ - + 1^1 V) + 2|M| 2 ^ 

+ + 2wf7 s/)l^| 2 " 2 K " (T a )ReT4 y (CT a + «)) + 

+ 5^Ao^(o2 + 4t«7- x )+w > (8) 

x 

where & = to — Xv 4 V. 

Now we easily derive the tree level estimates: 

M s = y/2\n\ = 2V4+(2A-l)/ 7 ; 



M fo v 7(4 7 + 2A-l) 

M w = cos9 w M z 

M H /M W = vW/7 2 ; 

A = 7T 



2A/3 



, — [80GeVl; (9) 
y 7(47 + 2A- 1) 1 J ' w 

The fine structure constant is given by a — j^^^ and does not depend on A and 7. From ([9]) we learn that at 
the tree level LCP on the phase diagram corresponds to fixed f3 = — n4r^fg — 7 ~ 10 and 77 = Mh/Mw, and is given 



by the equation A( 7 ) = |^7 2 - 

The important case is A = 00, where the tree level estimates give 




(10) 

In the SU(2) gauge Higgs model for the small values of A << 0.1 the tree level estimate for Mh/Mw gives values 
that differ from the renormalized ratio by about 20% ■ The tree level estimate for the ultraviolet cutoff is about 1 
TeV at A = oo,7 = 1, /3 = 15 that is not far from the numerical result given in 0. In the SU{2) Gauge Higgs model 
at A = 00 the critical 7 C = 0.63 for f3 = 8 (20| . At this point the tree level estimate gives A = 0.9 Tev while the direct 
measurements give A e [0.8; 1.5] Tev for values of 7 G [0.64; 0.95] [2(j. The investigations of the SU(2) Gauge Higgs 
model showed that a consideration of finite A does not change much the estimate for the gauge boson mass. However, 
at finite A and values of 7 close to the phase-transition point the tree level formula does not work at all. 

The tree level estimate for the critical 7 is -f c = (1 — 2A)/4. At small A this formula gives values that are close to 
the ones obtained by the numerical simulations [l8l - l20| . In particular, j c — > 0.25 (n c — > 0.125) at A << 1. However, 
this formula clearly does not work for A > 1/2. From [TsT - Ebl f30lj we know that the critical coupling in the SU{2) 
Gauge Higgs model is about 2 — 4 times smaller for A = than for A = 00. 

Tree level estimate predicts that there is the second order phase transition. This means that according to the tree 
level estimate the value of the cutoff at the transition point is infinite. Our numerical simulations, however, show that 
the cutoff remains finite and the transition is, most likely, a crossover at the considered values of 9yy, A and j3. 
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IV. NAMBU MONOPOLES 



In this section we remind the reader what is called Nambu monopole [2lj j . First let us define the continuum 
Electroweak fields as they appear in the Weinberg-Salam model. The continuum scalar doublet is denoted as $. The 
Z-boson field Z M and electromagnetic field A^ M are defined as 

= — — $+^'$ - B», 

,4£ M = 2^ + 2sin 2 %-Z", (11) 

where and B M are the corresponding SU(2) and U(l) gauge fields of the Standard Model. 
After fixing the unitary gauge $2 = const., $1 = we have 

2 92 9i 2 

A ^ B^ 

A£ M = e [— sin(9 w + cos0 w ] = & , (12) 

92 gi 



where Af = j^TrAa 3 , S. = 2Bjg x , j-, A - conventional Standard Model fields, and g z = + 92- 

Nambu monopoles are defined as the endpoints of the Z-string [2l| . The Z-string is the classical field configuration 
that represents the object, which is characterized by the magnetic flux extracted from the Z-boson field. Namely, for 
a small contour C winding around the Z - string one should have 

Z^dx** ~ 2vr; f A^dx** ~ 0; / - 27rsin 2 9 W . (13) 



JC 

The string terminates at the position of the Nambu monopole. The hypercharge flux is supposed to be conserved 
at that point. Therefore, a Nambu monopole carries electromagnetic flux 47rsin 2 9w- The size of Nambu monopoles 
was estimate d f2lll to be of the order of the inverse Higgs mass, while its mass should be of the order of a few TeV. 
According to [2l| Nambu monopoles may appear only in the form of a bound state of a monopole-antimonopole pair. 
In lattice theory the following variables are considered as creating the Z boson: 

Z xy = Z% = -sin [Arg^+E^e*^)]. (14) 

and: 

Z' xy = Z£ = - Anw .I-;/' ( "'-!.„) . (15) 

The classical solution corresponding to a Z-string should be formed around the 2-dimensional topological defect 
which is represented by the integer-valued field defined on the dual lattice S = A. ([dZ'] mo d2ir — dZ'). (Here we used 

the notations of differential forms on the lattice. For a definition of those notations see, for examplej 22[. Lattice 

field Z' is defined in Eg. (|15l) .) Therefore, S can be treated as the worldsheet of a quantum Z-string [23|. Then, the 
worldlines of quantum Nambu monopoles appear as the boundary of the Z-string worldsheet: jz = 8H. 

For historical reasons in lattice simulations we fix unitary gauge $2 = 0; $1 6 K; $1 > (instead of the usual 
$1 = 0; $2 G and the lattice Electroweak theory becomes a lattice U(l) gauge theory with the U(l) gauge field 

A*v = A x = i z ' + 2d *v] mod 27r . ( 16 ) 

(The usual lattice Electromagnetic field is related to A as A-^m = A — Z' + 2 sin 2 6wZ'.) One may try to extract 
monopole trajectories directly from A. The monopole current is given by 

OA = d{[dA]uvoA2Tr) (17) 

Both jz, and ja carry magnetic charges. That's why it is important to find the correspondence between them. 
In continuum notations we have 

A» = Z I * + 2B», (18) 

where B is the hypercharge field. Its strength is divergenceless. As a result in continuum theory the net Z flux 
emanating from the center of the monopole is equal to the net A flux. (Both A and Z are undefined inside the 
monopole.) This means that in the continuum limit the position of the Nambu monopole must coincide with the 
position of the antimonopole extracted from the field A. Therefore, one can consider Eq. (|17[) as another definition 
of a quantum Nambu monopole [3] . Actually, in our numerical simulations we use the definition of Eq. (1171) . 
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FIG. 1: The phase diagram of the model in the (7, A)-plane at ft = 12. The dashed line is the tree - level estimate for the 
line of constant physics correspondent to bare M% = 270 Gev. The continuous line is the line of phase transition between the 
physical Higgs phase and the unphysical symmetric phase (statistical errors for the values of 7 at each A on this line are about 
0.005). 





FIG. 2: Susceptibility x = (H 2 ) — (H) 2 (for H x = Y2, y Z 2 y ) as a function of 7 at A = 0.001 and ft = 12. Circles correspond to 



the lattice 12 3 x 16. Crosses correspond to the lattice 8 3 x 16. 



PHASE DIAGRAM 



dimensional (/3, 7, A) phase diagram the transition 

A _ 1 Ml 



w<m9- = const ; 



In our lattice study we fix bare 6w — 7r/6. Then in the three 

surfaces are two - dimensional. The lines of constant physics on the tree level are the lines (— ^ — — —r 

j3 = — = const) . We suppose that in the vicinity of the transition the deviation of the lines of constant physics from 
the tree level estimate may be significant. However, qualitatively their behavior is the same. Namely, the cutoff is 
increased along the line of constant physics when 7 is decreased and the maximal value of the cutoff is achieved at the 
transition point. Nambu monopole density in lattice units is also increased when the ultraviolet cutoff is increased. 

At (3 — 12 (corresponds to bare a ~ 1/150) the phase diagram is represented on Fig. Q] This diagram is obtained, 
mainly, using the lattice 8 3 x 16. Some regions (A = 0.009,0.0025,0.001), however, were checked using larger lattices. 
According to our data there is no dependence of the diagram on the lattice size. The physical Higgs phase is situated 
right to the transition line. The position of the transition 7c (A) is localized here at the point where the susceptibility 
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FIG. 3: The effective constraint potential at A = 0.009 and f3 = 12. Black squares correspond to -y c = 0.273. Empty squares 
correspond to 7 = 0.29. Triangles correspond to 7 = 0.279. The error bars are about of the same size as the symbols used. 



extracted from the Higgs field creation operator achieves its maximum. We use the susceptibility 

X =(H 2 )-(Hf (19) 

extracted from H = Z% y (see, for example, Fig. [5]). We observe no difference between the values of the suscepti- 
bility calculated using the lattices of different sizes. This indicates that the transition at 7 C is a crossover. Indeed we 
find that gauge boson masses do not vanish in a certain vicinity of j c even within the symmetric phase. In the next 
section we shall see that within the statistical errors j c coincides with the value of 7, where the scalar field condensate 
disappears. Actually, there also exist two other crucial points: 7 c o(A) < 7 C (A) < 7 C 2(A) (say, at A = 0.001 we have 
7 c o = 0.252 ± 0.001, 7 C = 0.256 ± 0.001, 7 c2 = 0.258 ± 0.001, see the next sections for the details). 7 c2 denotes the 
boundary of the fluctuational region. At j c0 the extrapolation of the dependence of lattice Z - boson mass Mz(j) 
on 7 indicates that Mz{^ c o) may vanish. In the symmetric phase the perturbation theory predicts vanishing of the 
gauge boson masses. Therefore, supposition that Mz vanishes at a certain point is very natural. The perturbation 
theory also predicts that the mass parameter present in the effective action for the scalar field vanishes at the point 
of the transition between Higgs phase and the symmetric phase. Our analysis shows that at the point, where the 
scalar field condensate disappears lattice M# does not vanish. However, it may vanish, in principle, at some other 
point. If both Mz and My vanish simultaneously at j c q, at this point the model becomes scale invariant and formal 
continuum limit of the lattice model can be achieved at "f C Q. This point may then appear as the point of the second 
order phase transition. Near "/ c q the fluctuations of the gauge boson correlator are strong and at the present moment 
we do not make definite conclusions on the behavior of the system at 7 c o- However, the calculated susceptibilities do 
not have peaks at this point that is an indirect indication that the real second order phase transition cannot appear 
at 7 c o. It is worth mentioning that within the region (7 c o,7c) the scalar field is not condensed. That's why we guess 
this region has nothing to do with real continuum physics. 

We investigated carefully the region 7 > j c for A = 0.001, 0.0025, 0.009. We observe that for 7 C < 7 < 7c2 Nambu 
monopoles dominate vacuum and the usual perturbation theory cannot be applied. For this reason, most likely, the 
interval (7 C ,7 C 2) also has no connection with the conventional continuum Electroweak theory. At the same time for 
7 >> 7c2 the behavior of the system is close to what one would expect basing on the usual perturbative continuum 
Weinberg - Salam model. It is worth mentioning that the value of the renormalized Higgs boson mass does not deviate 
significantly from its bare value near the transition point j c . For example, for A around 0.009 and 7 = 0.274 bare 
value of the Higgs mass is around 270 Gev while the observed renormalized value is 300 ± 70 Gev. 



VI. EFFECTIVE CONSTRAINT POTENTIAL 



We have calculated the constraint effective potential for |$| using the histogram method. The calculations have been 
performed on the lattice 8 3 x 16. The probability h(<fi) to find the value of |$| within the interval [<f> — 0.05; <p + 0.05) 
has been calculated for (j> = 0.05 + N * 0.1, N = 0, 1,2, ... This probability is related to the effective potential as 
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FIG. 4: 4> m as a function of 7 at A = 0.001 and f3 — 12. Circles correspond to lattice 8 3 x 16. Crosses correspond to lattice 
12 3 x 16. 
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FIG. 5: H (points) vs. Hfi uct (stars) as a function of 7 at A = 0.009 and /3 = 12. Statistical errors for Hfi uct are about of the 
same size as the symbols used. 



h(4>) = 4> 3 e That's why we extract the potential from h(<j>) as 

V(4>) = -logh(<p) + 3log<p (20) 

(See Fig. O) It is worth mentioning that /i(0.05) is calculated as the probability to find the value of |$| within the 
interval [0;0.1]. Within this interval \og(j> is ill defined. That's why we exclude the point <f> = 0.05 from our data. 
Instead we calculate V(0) using the extrapolation of the data at 0.15 < 4> < 2.0. The extrapolation is performed using 
the polynomial fit with the powers of <f> up to the third (average deviation of the fit from the data is around 1 per 
cent). Next, we introduce the useful quantity H = V(0) — V(<p m ), which is called the potential barrier hight (here 
4> m is the point, where V achieves its minimum). 

As an example we represent on Fig. 0]the values of </> m for A = 0.001, f3 = 12. On Fig. [5]we represent the values of 
H for A = 0.009, j3 = 12. One can see that the values of <t> m and H increase when 7 is increased. The maximum of 
the susceptibility constructed of the Higgs field creation operator H x = ^ Z% (see, for example, Fig. \2§ coincides 
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FIG. 6: Mean value of \(f>\ as a function of 7 at A = 0.0025 and /3 = 12. (Lattice 8 3 x 16.) 



with the point, where 4> m vanishes within the statistical errors. We localize the position of the transition points at the 
points where <f> m vanishes: 7 C = 0.274 ± 0.001 at A = 0.009; j c = 0.26 ± 0.001 at A = 0.0025; and 7 C = 0.256 ± 0.001 
at A = 0.001. 

The maximum of the scalar field fluctuation (see, for example, Fig. [7]) is shifted to larger values of 7 than the 
transition point. Again we do not observe any difference in 8(f) for the considered lattice sizes. This also indicates 
that the transition at these values of A is a crossover. 

It is important to understand which value of barrier hight can be considered as small and which value can be 
considered as large. Our suggestion is to compare H — V(0) — V(<j> m ) with iffl uct = V((f> m + 8(f)) — V((f> m ), where 
8(f> is the fluctuation of |$|. From Fig. [5] it is clear that there exists the value of 7 (we denote it 7 c2 ) such that at 
7 C < 7 < 7 C 2 the barrier hight H is of the order of -fffl uct while for 7 C 2 << 7 the barrier hight is essentially larger 
than iffluct- The rough estimate for this pseudocritical value is 7 C 2 ~ 0.278 at A = 0.009. 

The fluctuations of |$| are around 8(f) ~ 0.6 for all considered values of 7 at A = 0.009,0.0025,0.001, /? = 12. It 
follows from our data (see also Fig. [6]) that c/> m , (\(f)\) » S(f> at 7 c2 << 7 while </> m , (\(f>\) ~ S(f> at 7 c2 > 7. Basing on 
these observations we expect that in the region 7 C 2 << 7 the usual perturbation expansion around trivial vacuum of 
spontaneously broken theory can be applied to the lattice Weinberg - Salam model while in the FR j c < 7 < 7 c2 it 
cannot be applied. In the same way we define the pseudocritical value 7 c2 at A = 0.001, 0.0025. Namely, 7 c2 ~ 0.278 
for A = 0.009; ~ 0.262 for A = 0.0025; ~ 0.258 for A = 0.001. 



VII. THE RENORMALIZED COUPLING 



In order to calculate the renormalized fine structure constant ccr — e 2 /47r (where e is the electric charge) we use 
the potential for infinitely heavy external fermions. 

We consider Wilson loops for the right-handed external leptons: 

< pt (0HRen ( ^ )ei e 2 ^). (21) 

Here I denotes a closed contour on the lattice. We consider the following quantity constructed from the rectangular 
Wilson loop of size r x t: 

V(r)=log lim , W( 7*> . (22) 
v ' 8 t^oo W(r x (i + 1)) v ' 

Due to exchange by virtual photons at large enough distances we expect the appearance of the Coulomb interaction 

V(r) = -— + const. (23) 
r 
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FIG. 7: Fluctuation 8<j> as a function of 7 at A = 0.0025 and f3 = 12. Circles correspond to the lattice 8 3 x 16. Crosses 
correspond to the lattice 12 3 x 16. Triangles correspond to the lattice 16 4 . Transition point is "f c — 0.261 ± 0.001; it is clear 
that the maximum of 5<f> is shifted to larger values of 7. 




FIG. 8: The potential (T = 8) for the right - handed leptons vs. 1/R at 7 = 0.262, A = 0.0025, and = 12 (lattice 16 4 ). 



It should be mentioned here, that in order to extract the renormalized value of a one may apply to V the fit obtained 
using the Coulomb interaction in momentum space. The lattice Fourier transform then gives 

V(r) = —aftU(r) + const, 

7r ^ e ip3T 
U{r) = — Y^^- 'L— (24) 

p^o sm Pi/2 + sin p 2 /2 + sin P3/2 

Here N is the lattice size, pi — =r-ki, fcj = 0, L—l. On large enough lattices at r << L both definitions approach 
each other. On the lattices we use the values of the renormalized Or extracted from (|2Uf and (|2~4")l are essentially 
different from each other. Any of the two ways, (|23p or (|24p. may be considered as the definition of the renormalized 
a on the finite lattice. And there is no particular reason to prefer the potential defined using the lattice Fourier 
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FIG. 9: The inverse renormalized fine structure constant as a function of 7 at A = 0.0025 and /3 = 12. It slowly approaches the 
tree level estimate ~ 150 when 7 and the lattice size are increased. Circles correspond to lattice 12 3 x 16. Crosses correspond 
to lattice 8 3 x 16. Triangles correspond to lattice 16 4 . 



transform of the Coulomb law in momentum space. Actually, our study shows that the single 1/r fit approximates V 
much better. Moreover, the values of renormalized a calculated using this fit are essentially closer to the tree level 
estimate than that of calculated using the fit ([24)) . 

In practise instead of (|22p we use the potential that depends on additional parameter T: 

For example, on the lattice 16 4 the values T — 4, 5, 6, 7, 8 are used; on the lattice 12 3 x 16 the values T = 4,5,6 are 
used; on the lattice 8 3 x 16 the value T = 4 is used. As a result a R = or(T) may depend both on the lattice size and 
on T. The dependence on T was missed in (where for lattices 12 3 x 16, 16 4 we used T = 5, while for the lattice 
8 3 x 16 we used T = 4). 

On Fig. [5] we represent as an example the dependence of the potential for T = 8 on 1 /R. As it was already 
mentioned (|23|) approximates the potential much better than ([24]) . Therefore we used the fit ([23]) to extract or. This 
should be compared with the results of 20], where for similar reasons the single e _Air /r fit (instead of the lattice 
Yukawa fit) was used in order to determine the renormalized coupling constant in the SU (2) Gauge Higgs model. 

Due to the dependence of cur (T) on T there is the essential uncertainty in definition of a R related to finite volume 
effects. For example, at 7 = 0.29, A = 0.009, and (3 = 12 the value of an calculated on the lattice 16 4 varies between 
Q!fl(4) ~ 1/(93±1) and a R (8) ~ 1/(108±2) (at the same time on the lattice 8 3 x 16 the value is a R {4) = 1/(100±1)). 
At 7 = 0.274, A = 0.009, and /3 = 12 the value of or calculated on the lattice 20 3 x 24 varies between or (4) ~ 1/(98±1) 
and or(IO) = 1/(106 ± 1) (at the same time on the lattice 8 3 x 16 the value is c*r(4) = 1/(99 ± 1)). Below for the 
lattice 8 3 x 16 we use T = 4, for the lattice 12 3 x 16 we use T = 6, for the lattice 16 4 we use T = 8. Therefore, 
the dependence on T is absorbed into the dependence on the lattice size. As an example, on Fig. [5] we represent the 
renormalized fine structure constant (calculated using the fit ([2"5]0 at A = 0.0025, /? = 12. The calculated values are 
to be compared with bare constant ao — 1/(4^/3) ~ 1/150 at j3 = 12. One can see, that for 7 >> 7 C 2 the tree level 
estimate is approached slowly while within the FR the renormalized a differs essentially from the tree level estimate. 
This is in correspondence with our supposition that the perturbation theory cannot be valid within the FR while it 
works well far from the FR. The dependence of a R on the lattice size is clear: for the larger lattices a R approaches its 
tree level estimate faster than for the smaller ones. Unfortunately, due to the difficulties in simulation of the system 
at large 7 we cannot observe this pattern in detail. At the present moment the value of an most close to the tree 
level estimate is obtained on the lattice 12 3 x 16 and is about 1/140 (at A = 0.0025, 0.001; (3 — 12; 7 ~ 1). 
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FIG. 10: Z - boson mass in lattice units at A = 0.009 and /3 = 12 as a function of 7. Black triangles correspond to lattice 
12 3 x 16. Crosses correspond to lattice 8 3 x 16. Circles correspond to lattice 16 4 . Square corresponds to lattice 20 3 x 24. (The 
error bars for lattices 16 3 x 16 and 20 3 x 24 are about of the same size as the symbols used.) 



VIII. MASSES AND THE LATTICE SPACING 



After fixing the unitary gauge $1 € R, $2 = 0, $1 > the following variables are considered as creating a Z boson 
and a W boson, respectively: 



Z xy = Z% = -sin[Arg(£^e"-»)] 

W xy = W£ = Ul 2 y e- ie *v. (26) 
Here, /i represents the direction (xy). The electromagnetic U(l) symmetry remains: 

U X y 9 x Uxy9yi 
@xy ^ @xy 

a y /2 + a x /2, (27) 

where g x — diag(e Ma =/ 2 , e~ la =°^ 2 ). There exists a f7 (1) lattice gauge field, which is defined as 

A xy = A% = [-Argt/ii + 6 xy ] mod 2tt (28) 

that transforms as A xy — > A xy — a y + a x . The field W transforms as W xy — > W xy e~ la:c . 

The W boson field is charged with respect to the U(l) symmetry. Therefore we fix the lattice Landau gauge in 
order to investigate the W boson propagator. The lattice Landau gauge is fixed via minimizing (with respect to the 
U(l) gauge transformations) the following functional: 

F = - cos(A xy )). (29) 

xy 

Then we extract the mass of the W boson from the correlator 

TV? EE W x( w v) ] ) ~ e- Mwlx °- yo1 + e- Mw{L -\ x °- y °V (30) 

x,y p. 

Here the summation - is over the three "space" components of the four - vectors x and y while Xo, yo denote their 
"time" components. N is the lattice length in "space" direction. L is the lattice length in the "time" direction. 
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FIG. 11: Z - boson mass in lattice units at A = 0.001 and /3 — 12. Circles correspond to lattice 8 3 x 16. Crosses correspond to 
lattice 12 3 x 16. 




In order to evaluate the masses of the Z-boson and the Higgs boson we use the correlators: 

_L Z x Z v) ~ e- Mz \ x °- y °\ + e -^(i- |xo- OT |) ( 31 ) 

x,y n 

and 

^2((H x H y ) - (H) 2 ) ~ e - M H|xo- W | + e -M H (L-\x -y \)^ ( 32 ) 
x,y 

In lattice calculations we used two different operators that create Higgs bosons: H x = |$| and H x = J2 y %xy ^ n 
both cases H x is defined at the site x, the sum J2 y i s over lts neighboring sites y. 

The physical scale is given in our lattice theory by the value of the Z-boson mass M v z ys ~ 91 GeV. Therefore the 
lattice spacing is evaluated to be a ~ [91GeV] _1 M^, where Mz is the Z boson mass in lattice units. The similar 
calculations have been performed in Q for A = oo. It has been found that the W - boson mass contains an artificial 
dependence on the lattice size. We suppose, that this dependence is due to the photon cloud surrounding the W - 
boson. The energy of this cloud is related to the renormalization of the fine structure constant. Therefore the Z - 
boson mass was used in order to fix the scale. 

Our data show that A = - = (it x 91 GeV)/Mz is increased slowly with the decrease of 7 at any fixed A. We 
investigated carefully the vicinity of the transition point at fixed A = 0.001,0.0025,0.009 and f3 = 12. It has been 
found that at the transition point the value of A is equal to 1.4 ± 0.2 TeV for A = 0.009, 0.0025, 0.001. Check of the 
dependence on the lattice size (8 3 x 16, 12 3 x 16, 16 4 , 20 3 x 24 at A = 0.009; 8 3 x 16, 12 3 x 16, 16 4 at A = 0.0025; 8 3 x 16, 
12 3 x 16 at A = 0.001) does not show an essential dependence of this value on the lattice size. This is illustrated 
by Fig. [TOl FigfTTI and Fig. [12] From these figures it also follows that at the value of 7 equal to 7 C 2(~ 0.278 for 
A = 0.009; - 0.262 for A = 0.0025; - 0.258 for A = 0.001) the calculated value of the cutoff is about 1 TeV. 

It is worth mentioning that the linear fit applied (in some vicinity of j c ) to the dependence of Mz on 7 predicts 
vanishing of Mz(j) at 7 equal to j c q < j c . Within the statistical errors j c o = 0.253 ± 0.001 for A = 0.001, 
7 c0 = 0.253 ± 0.001 for A = 0.0025, 7 c0 = 0.254 ± 0.001 for A = 0.009. We perform direct calculations within the 
region (7 c o,7c) at A = 0.001,0.0025. These calculations show that the fluctuations of the correlator (f3"Tj) are increased 
(compared with the values of the correlator) fast when 7 is decreased. Already for 7 = 0.255 at A = 0.0025 (<y c — 0.26) 
and for 7 = 0.254 at A = 0.001 (j c = 0.258) the values of the correlator at |xo — yo\ > are smaller than the statistical 
errors. Most likely, at 7 < j c0 it is necessary to apply another gauge (like in pure SU(2) x {7(1) gauge model) in 
order to calculate gauge boson propagators. At the present moment we do not estimate the scalar particle mass at 7 c o 
because of the lack of statistics. The behavior of the other quantities is smooth at 7 ~ 7 c o, no maximum of 6(j> or other 
susceptibilities is observed there (see, for example, Fig. [3]). Basing on our data it is natural to suppose that lattice 
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FIG. 12: Z - boson mass in lattice units at A = 0.0025 and j3 = 12. Circles correspond to lattice 12 3 x 16. Crosses correspond 
to lattice 8 3 x 16. Triangles correspond to lattice 16 4 . 



gauge boson mass may vanish at 7 ~ j c q although we do not observe the correspondent pattern in details because 
of the strong fluctuations of correlator (|31|) near 7 c o. As it was mentioned above the transition for the considered 
values of couplings is, most likely, a crossover. There are 3 exceptional points: 7 c o, where lattice value of Mz may 
vanish, 7 C , where scalar field condensate disappears, and 7 C 2 that denotes the boundary of the fluctuational region. 
This situation is typical for the crossovers: different quantities change their behavior at different points on the phase 
diagram. At the present moment we do not exclude that the second order phase transition may take place at 7 c o- 
This would happen if both mass parameters (Z boson mass and scalar particle mass) vanish simultaneously at this 
point. The careful investigation of the vicinity of 7 c o is to be the subject of a further research. 

In the Higgs channel the situation is more difficult. Due to the lack of statistics we cannot estimate the masses 
in this channel using the correlators (J32J) at all considered values of coupling constants. Moreover, at several points, 
where we have estimated the rcnormalized Higgs boson mass the statistical errors are much larger than that of for the 
Z - boson mass. At the present moment we can represent the data at four points on the lattice 8 3 x 16: (7 = 0.274, 



A = 0.009, (3 = 12), (7 = 0.290, A = 0.009, /3 = 12), (7 = 0.261, A = 0.0025, (3 = 12), and (7 = 0.257, A = 0.001, 



The first point roughly corresponds to the position of the transition at A = 0.009, f3 = 12 while the second point 
is situated deep within the Higgs phase. These two points correspond to bare Higgs mass around 270 Gev. At the 
point (7 = 0.274, A = 0.009, j3 = 12) we have collected enough statistics to calculate correlator (1521 up to the "time" 
separation \xq — 2/0 1 =4. The value 7 = 0.274 corresponds roughly to the position of the phase transition. We 
estimate at this point Mh = 300 ± 40 Gev. At the point (7 = 0.29, A = 0.009, j3 = 12) we calculate the correlator 
with reasonable accuracy up to \xq — yo| =3. At this point Mh = 265 ± 70 Gev. 

For A = 0.001, 0.0025 we calculate the Higgs boson mass close to the transition points. Similar to the case A = 0.009 
we do not observe here essential deviation from the tree level estimates. Namely, for A = 0.001,7 = 0.257 we have 
Mh = 90 ± 20 GeV (tree level value is M^ ~ 100 GeV). In this point we have collected enough statistics to calculate 
correlator (J32J) up to the "time" separation \x - yo\ = 8. For A = 0.0025,7 = 0.261 we have M H = 170 ± 30 GeV 
(tree level value is M H ~ 150 GeV). In this point we have collected enough statistics to calculate correlator (|3"2"]) up to 
the "time" separation |jcq — yo\ =4. It is worth mentioning that in order to calculate Z - boson mass we fit correlator 
HE) for 8 > \x - y \ > 1. 



The worldlines of the quantum Nambu monopoles can be extracted from the field configurations according to Eq. 
dl~7|) . The monopole density is defined as p = ( ^ lin ^^''' nk l \ 5 where V L is the lattice volume. 



On Fig[T31 Fig. Q31 and Fig. Ql)]we represent Nambu monopole density as a function of 7 at A = 0.009, 0.0025, 0.001, 



(3 = 12). 



IX. NAMBU MONOPOLE DENSITY 
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FIG. 13: Nambu monopole density as a function of 7 at A = 0.009 and (3 — 12. (Lattice 8 3 x 16.) 



(5 = 12. The value of monopole density at 7 C is around 0.1. 

According to the classical picture the Nambu monopole size is of the order of Mjj . Therefore, for example, for 
a -1 <~ 430 Gev and M H <~ 300,150,100 Gev the expected size of the monopole is about a lattice spacing. The 
monopole density around 0.1 means that among 10 sites there exist 4 sites that are occupied by the monopole. 
Average distance between the two monopoles is, therefore, less than 1 lattice spacing and it is not possible at all to 
speak of the given configurations as of representing the physical Nambu monopole. 

At 7 = 7 C 2 the Nambu monopole density is of the order of 0.01. This means that among about 25 sites there exists 
one site that is occupied by the monopole. Average distance between the two monopoles is, therefore, between one 
and two lattice spacings. We see that at this value of 7 the average distance between Nambu monopoles is of the 
order of their size. 

We summarize the above observations as follows. Within the fluctuational region the configurations under consid- 
eration do not represent single Nambu monopoles. Instead these configurations can be considered as the collection 
of monopole - like objects that is so dense that the average distance between the objects is of the order of their size. 
On the other hand, at 7 >> 7 C 2 the considered configurations do represent single Nambu monopoles and the average 
distance between them is much larger than their size. In other words out of the FR vacuum can be treated as a gas 
of Nambu monopoles while within the FR vacuum can be treated as a liquid composed of monopole - like objects. 

It is worth mentioning that somewhere inside the Z string connecting the classical Nambu monopoles the Higgs 
field is zero: |$| = 0. This means that the Z string with the Nambu monopoles at its ends can be considered as an 
embryo of the symmetric phase within the Higgs phase. We observe that the density of these embryos is increased 
when the phase transition is approached. Within the fluctuational region the two phases are mixed, which is related 
to the large value of Nambu monopole density. 

That's why we come to the conclusion that vacuum of lattice Weinberg - Salam model within the FR has nothing to 
do with the continuum perturbation theory. This means that the usual perturbation expansion around trivial vacuum 
(gauge field equal to zero, the scalar field equal to (4> m , 0) T ) cannot be valid within the FR. This might explain why we 
do not observe in our numerical simulations the large values of A predicted by the conventional perturbation theory. 

X. CONCLUSIONS 

In the present paper we demonstrate that while approaching continuum physics in lattice Weinberg - Salam model 
one encounters the nonperturbative effects. Namely, the continuum physics is to be approached in the vicinity of the 
transition between the physical Higgs phase and the symmetric phase of the model (in the symmetric phase the scalar 
field is not condensed). The ultraviolet cutoff is increased when the transition point is approached along the line of 
constant physics. There exists the fluctuational region (FR) on the phase diagram of the lattice Weinberg - Salam 
model. This region is situated in the vicinity of the transition between the Higgs phase and the symmetric phase 
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FIG. 14: Nambu monopole density as a function of 7 at A = 0.0025 and f3 = 12. Circles correspond to lattice 12 x 16. Crosses 
correspond to lattice 8 3 x 16. Triangles correspond to lattice 16 4 . 
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FIG. 15: Nambu monopole density as a function of 7 at A = 0.001 and j3 = 12. Circles correspond to lattice 12 3 x 16. Crosses 
correspond to lattice 8 3 x 16. 



(where scalar field is not condensed). According to our data this transition is, most likely, a crossover. We localize its 
position at the point 7 C (A, /3, 6w), where the scalar field condensate disappears. We calculate the effective constraint 
potential V{<j>) for the Higgs field. It has a minimum at the nonzero value 4> m in the physical Higgs phase. At the 
considered values of A, (3, 9w for 7 between 7 C and 7 c2 (7 C 2 is in the Higgs phase) the fluctuations of the scalar field 
become of the order of <fi m . Moreover, the "barrier hight" H = V(0) — V((p m ) is of the order of V{<j> m + 8<j>) — V(<f> m ), 
where 5(f> is the fluctuation of |$|. Therefore, we refer to this region as to FR. 

The scalar field must be equal to zero somewhere within the classical Nambu monopole. That's why this object 
can be considered as an embryo of the unphysical symmetric phase within the physical Higgs phase of the model. We 
investigate properties of the quantum Nambu monopoles. Within the FR they are so dense that the average distance 
between them becomes of the order of their size. This means that the two phases are mixed within the FR. All these 
results show that the vacuum of lattice Weinberg - Salam model in the FR is essentially different from the trivial 
vacuum used in the conventional perturbation theory. As a result the use of the perturbation theory in this region is 
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limited. 

Our numerical results show that at M H around 270, 150, 100 GeV and the bare fine structure constant around 
1/150 the maximal value of the cutoff admitted out of the FR for the considered lattice sizes cannot exceed the value 
around 1 Tev. Within the FR the larger values of the cutoff can be achieved in principle. The maximum for the value 
of the cutoff A c within the Higgs "phase" of the lattice model is achieved at the point of the transition to the region 
of the phase diagram, where the scalar field is not condensed. Our estimate for this value is A c = 1.4 ± 0.2 Tev for 
the considered lattice sizes. Far from the fluctuational region the behavior of the lattice model in general is close to 
what we expect basing on the continuous perturbation theory. As it was already mentioned at the considered values 
of couplings the transition is, most likely, a crossover. This follows from the observation that various quantities (Z 
boson mass, the fluctuation of the scalar field etc) do not depend on the lattice size at the transition point. Within 
the symmetric " phase" of the lattice model (where the scalar field is not condensed) in some vicinity of the transition 
between this phase and the Higgs phase (where the scalar field is condensed) the lattice gauge boson masses do not 
vanish. The statistical error for Mz is increased fast when 7 is decreased starting from the pseudocritical value j c . At 
7 < 7c0 < 7c (within the symmetric phase) the values of the Z - boson correlator ([3Tj) are smaller than the statistical 
errors. Therefore, our procedure cannot give the values of gauge boson masses in this region. Most likely, here the 
other gauge is to be applied in order to calculate gauge boson propagators (we used in our simulations the Unitary 
gauge). It is worth mentioning that the perturbation theory predicts zero gauge boson masses within the symmetric 
phase. Most likely, this prediction is failed within the interval (7 C 0i7c) due to nonperturbative effects. 

An important question is how to treat finite volume effects that are present in all observables that contain long - 
ranged Electromagnetic Coulomb interactions. In particular, we see that these effects are strong in renormalized fine 
structure constant (about 10% when the lattice size varies from 8 3 x 16 to 16 4 ) and in the mass of electrically charged 
W - boson. On the other hand all observables related to SU (2) constituent of the model do not possess essential 
dependence on the lattice size. In particular, Z - boson mass Mz (and the cutoff A), density pNamtm of Nambu 
monopoles 39], fluctuation of the scalar field 5cj) as well as the position of the transition between the "phases" of the 
lattice model practically do not depend on the lattice size. Our point of view is that the influence of long - ranged 
Electromagnetic interactions on these observables is negligible compared to their tree - level and nonperturbative 
constituents. Actually, Electromagnetic interactions can be taken into account perturbatively, with the renormalized 
a r~j 1/100 as the parameter of the perturbation expansion. This was the reason why in the previous numerical studies 
of SU{2) Gauge - Higgs model the U(l) constituent of Weinberg - Salam model was completely disregarded 0413 ■ 
To summarize, we suppose that in spite of the presence of finite volume effects in fine structure constant and W 
boson mass, the calculated values of Mz , A, pNambu, S(f> etc can be considered as free of these effects[40( (up to the 
perturbations suppressed by the factor a ~ 1/100). 

Basing on our data it is natural to suppose that lattice gauge boson mass may vanish at 7 ~ 7 c o although we 
do not observe the correspondent pattern in details because of the strong fluctuations of correlator (|31l) near 7 c o. If 
so, there exist 3 pseudocritical points: 7 c0 , where lattice value of Mz vanishes (at this point the cutoff calculated 
as A = (n x 91 GeV)/Mz tends to infinity), 7 C , where scalar field condensate disappears, and 7 c2 that denotes the 
boundary of the fluctuational region (at 7 ~ 7 C 2 the average distance between Nambu monopoles becomes of the 
order of their size). This situation is typical for the crossovers: different quantities change their behavior at different 
points on the phase diagram. There still exists the possibility that the point j c q corresponds to the second order 
phase transition (this may happen if, in addition, the scalar particle mass vanishes at 7 c o)- However, the absence 
of a peak in the scalar field fluctuation and in susceptibility (|19p at this point indicates that this is a crossover. 
Actually, this possibility is to be checked carefully but this is to be a subject of another work. There is an important 
question: what is the relation between the conventional Electroweak physics and the regions {^ c q,1c) and {^ c ,lc2)- 
Our expectation is that both these regions have nothing to do with real continuum physics. For the first region this is 
more or less obvious: there the scalar field is not condensed that contradicts with the usual spontaneous breakdown 
pattern. As for the second region, the situation is not so obvious. However, there the nonperturbative effects are 
strong and the Nambu monopoles dominate vacuum that seems to us unphysical. With all mentioned above we come 
to the conclusion that our data indicate the appearance of the maximal value of the cutoff in Electroweak theory that 
cannot exceed the value of the order of 1 TcV. This prediction is made basing on the numerical investigation of the 
lattice model on the finite lattices. However, as it was mentioned above, our main results do not depend on the lattice 
size. 



Acknowledgments 

This work was partly supported by RFBR grants 09-02-00338, 08-02-00661, by Grant for leading scientific schools 
679.2008.2. The numerical simulations have been performed using the facilities of Moscow Joint Supercomputer 



19 



Center. 



[1] Peter Arnold and Olivier Espinosa, Phys. Rev. D 47, 3546 (1993) 
Z. Fodor and A. Hebecker, Nucl. Phys. B 432, 127 (1994) 
W. Buchmuller, Z. Fodor, and A. Hebecker, Nucl. Phys. B 447, 317 (1995) 
B.L.G. Bakker, A.I. Veselov, M.A. Zubkov. J.Phys.G36, 075008 (2009) 
M.A. Zubkov, A.I. Veselov. JHEP 0812, 109 (2008) 
M.A. Zubkov. Phys.Lett.B 684, 141 (2010) 

K.Holland, Plenary talk presented at Lattice2004, Fermilab, June 21-26, 2004, |arXiv:hep-lat/0409112"1 
Zoltan Fodor, Kieran Holland, Julius Kuti, Daniel Nogradi, Chris Schroeder, The XXV International Symposium on 
Lattice Field Theory, July 30 - August 4 2007, Regensburg, Germany, PoS (LATTICE 2007), 056. larXiv:0710.3151l 
J. A. Casas, J.R. Espinosa, and I. Hidalgo, Nucl.Phys.B 777, 226 (2007) 

F. Csikor, Z. Fodor, J. Heitger Phys.Rev.Lett. 82, 21 (1999), Phys.Rev. D58, 094504 (1998), Nucl.Phys.Proc.Suppl. 63, 
569 (1998) 

F. Csikor, Z. Fodor, J. Heitger Phys.Lett. B 441, 354 (1998) 

F. Csikor, Z. Fodor, J. Hein, A. Jaster, I. Montvay Nucl.Phys. B 474, 421 (1996) 
Joachim Hein (DESY), Jochen Heitger, Phys.Lett. B 385, 242 (1996) 

F. Csikor, Z. Fodor, J. Hein, J. Heitger, A. Jaster, I. Montvay Nucl.Phys.Proc.Suppl. 53, 612 (1997) 
Z. Fodor, J. Hein, K. Jansen, A. Jaster, I. Montvay Nucl.Phys. B 439 (1995) 
F. Csikor, Z. Fodor, J. Hein, J. Heitger, Phys.Lett. B 357, 156 (1995) 

F. Csikor, Z. Fodor, J. Hein, K. Jansen, A. Jaster, I. Montvay Nucl.Phys.Proc.Suppl. 42, 569 (1995) 
F. Csikor, Z. Fodor, J. Hein, K. Jansen, A. Jaster, I. Montvay Phys.Lett. B 334, 405 (1994) 
Y. Aoki, F. Csikor, Z. Fodor, A. Ukawa Phys.Rev. D 60, 013001 (1999), Nucl.Phys.Proc.Suppl. 73, 656 (1999) 
Y. Aoki Phys.Rev. D 56, 3860 (1997) 

W.Langguth, I. Montvay, P.Weisz Nucl.Phys.B 277, 11 (1986) 
W. Langguth, I. Montvay, Z.Phys.C 36, 725 (1987) 
Anna Hasenfratz, Thomas Neuhaus, Nucl.Phys.B 297, 205 (1988) 
Y. Nambu, Nucl.Phys. B 130, 505 (1977); 

Ana Achucarro and Tanmay Vachaspati, Phys. Rept. 327, 347 (2000); Phys. Rept. 327, 427 (2000) 
M.I. Polikarpov, U.J. Wiese, and M.A. Zubkov, Phys. Lett. B 309, 133 (1993) 
M.N. Chernodub, JETP Lett. 66, 605 (1997) 

Bohdan Grzadkowski, Jose Wudka, Acta Phys. Polon. B 32, 3769 (2001) 

B.L.G. Bakker, A.I. Veselov, M.A. Zubkov. Yad.Fiz.68, 1045 (2005), Phys.Atom.Nucl.68, 1007 (2005) 
B.L.G. Bakker, A.I. Veselov, M.A. Zubkov. Phys.Lett.B 620, 156 (2005) 
B.L.G. Bakker, A.I. Veselov, M.A. Zubkov. Phys.Lett.B 642, 147 (2006) 

A.I. Veselov, B.L.G. Bakker, M.A. Zubkov, The XXV International Symposium on Lattice Field Theory, July 30 - August 
4 2007, Regensburg, Germany, PoS (LATTICE 2007), 337 (2007). larX"iv:0708. 28641 
I. Montvay, Nucl. Phys. B 269, 170 (1986) 
W.Langguth, I. Montvay, P.Weisz, Nucl.Phys.B 277, 11 (1986) 
R. Shrock, Phys. Lett. B 162, 165 (1985); Nucl. Phys. B 267, 301 (1986) 
I. Montvay, DESY preprint 86-143 (1986), DESY preprint 87-019 (1987) 
Bohdan Grzadkowski, Jacek Pliszka, Jose Wudka Phys.Rev. D 69, 033001 (2004) 
M.N. Chernodub, Phys.Rev.Lett. 95, 252002 (2005) 

According to the previous investigations of the SU (2) Gauge - Higgs model this upper bound cannot exceed 10Mw . 
One of the examples of such models is the Ginzburg - Landau theory of superconductivity. 

The meaning of the words "potential barrier" here is different from that of the one - dimensional quantum mechanics as 
here different minima of the potential form the three - dimensional sphere while in usual ID quantum mechanics with 
the similar potential there are two separated minima with the potential barrier between them. Nevertheless we feel it 
appropriate to use the chosen terminology as the value of the "potential barrier hight" measures the difference between 
the potentials with and without spontaneous symmetry breaking. 
[38] It has been shown in [3J that at the infinite value of the scalar self coupling A = oo moving along the line of constant 
physics we reach the point on the phase diagram where the monopole worldlines begin to percolate. This point was found 
to coincide roughly with the position of the transition between the physical Higgs phase and the unphysical symmetric 
phase of the lattice model. This transition is a crossover and the ultraviolet cutoff achieves its maximal value around 1.4 
Tev at the transition point. 

[39] Nambu monopoles in practise correspond to SU (2) variables as the monopole configurations extracted from the Hypercharge 
U(l) field disappear at realistic values of coupling constants. 

[40] The inverse seem to us incorrect: influence of nonperturbative effects on cxr is not suppressed by any small factor. We 
indeed observe that in the FR, where nonperturbative effects are large the renormalized a differs from its bare value by 
about 50% while far from the FR the difference is within 10% (for the lattice size 12 3 x 16). 



